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The NN scattering near inelastic threshold is sensitive to the long-range diagonal interaction in 
the produced isobar channel. Earlier models included meson exchange potentials in the NN sector 
and connecting that sector to isobar channels, as well as an .R-matrix description of short-range 
quark effects. Including the diagonal pion exchange contributions in the NA channels coupled to 
the very inelastic D2 and F3 NN channels substantially improves those phase parameters and the 
observables Acr^ and Ao~ T . The same improvement is made to the other 1=1 NN channels. In 
the 5q NN channel coupling to the A A channel is unusually important due to the L = 2 angular 
momentum barrier in the AA channel. The pion exchange transition potentials between isobar 
channels are included in this partial wave to obtain the correct equilibrium between one- and two- 
pion decay channels. Other improvements to earlier models have been made, in the specification of 
isobar channels and the inclusion of decay width effects in more channels. A comparison is made 
with all the pp data for Tl < 800 MeV, producing a very good fit over the whole energy range. 
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I. Introduction 

In describing the nucleon-nucleon (AW) interaction beyond the first 150 MeV of kinetic energy 
in the center of momentum system (nucleon beam energy Tl > 300 MeV), coupling to isobar 
channels introduces an important energy dependence to the real phase even when the imaginary 
phase remains small. Single-pion production begins at Tl = 280 MeV, but is relatively weak until 
the A production threshold, Tl = 630 MeV is approached. But the effect on the elastic scattering 
begins below the threshold of substantial inelasticity because of the effect of virtual intermediate 
states of higher mass. Consequently, a coupled channel system must be used to understand the 
NN reaction in the whole range < Tl < 800 MeV. Experience has shown 1 ' 2 ' 3 that for some NN 
partial waves the AA and AW* (1440) channels are important, in addition to the low threshold AA 
channel. 

The i?-matrix (or Boundary Condition) model 4 ' 5 incorporates long-range meson exchange ef- 
fects in the potential matrix outside the boundary radius tq. 1 The internal dynamics, including 
the effects of quark degrees of freedom, is described by a simple boundary condition at tq. The 
boundary condition is a meromorphic function of energy 4 whose poles are completely described by 
the properties of simple quark configurations. 2 ' 3 ' 6 ' 7 Only the constant term is not given by the 
internal model. 

Using a potential matrix and a boundary condition matrix the coupling to isobaric channels is 
easily incorporated. 1 With some computational burden, but little formal complication, the width of 
the isobars is included 8 by discretizing the Breit-Wigner mass distribution of the isobar, treating 
each mass as a separate channel. In previous work, 3, 8 this has provided a good quantitative descrip- 
tion of the phase parameters of most of the proton-proton (pp) partial waves. The most important 
exceptions to this quantitative agreement have been the detailed shapes of the inelasticity, rj, in 
the x T)<i and 3 F3 channels. 8 The predicted positions, strengths and widths of the structure in these 
channels are in agreement with experiment, indicating that the structures are due to the effect of 
the AA threshold and require no significant contribution from quark substructure (the latter is 
important at somewhat higher energies 2, 3 ' 6,7 ). In fact this potential and i?-matrix model has been 
the only one to adequately predict the strong onset of inelasticity in the 3 ^3 channel. However, 
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because of the dominance of these channels in the Aa L (pp) and Aa T (pp) structures, for 600 MeV 
< Tl < 800 MeV, 9 the detailed shaped of the inelasticity must be fairly precise to reproduce the 
data in the region. 

In the model, the 7r, 77, p and uj meson exchanges (and also two-pion exchange) are all included 
in deriving the potential within the AW sector (Fig. l(a)-(c)). 

Because of the indirect effect of the AW to isobar channel transition (Fig. l(d)-(e)) on the 
NN elastic interaction only the n meson exchange transition potentials are included. A two-pion 
exchange transition potential is only important at distances less than one-third of a pion Compton 
wavelength in the NN — > NN reactions, at which distance it is partly masked by the effect of 
the boundary condition. Nevertheless, a phenomenological two-pion range transition potential is 
sometimes included, which may also substitute for p meson exchange. 

The diagonal NA interaction is one stage further removed than the transition interaction from 
influencing the NN scattering, and was entirely represented by the boundary condition in previous 
work. However, when the energy dependence of the effect is strong, as it is near the NA threshold, 
then the short-range boundary condition is not an adequate substitute for the long-range diagonal 
potential. The long-range interaction shown in Fig. l(f)-(g) is more important at small momentum 
transfers. In this paper we include these diagonal potentials, improving many partial wave phases. 
In particular, the rj and 5 in the x D<i and 3 Fs channels are greatly improved for Tl > 500 MeV. 

In previous work 2 ' 3 there was some difficulty in obtaining a very good fit to the higher energy 

NN phases simultaneously with correctly predicting the experimentally very small two-pion 
production at Tl = 800 MeV. 10 The NA diagonal potential affects this situation positively, but 
the transition between NA channels, which produce one pion, and AA or AW* (1440) channels, 
which produce two pions, critically affects the pion production multiplicity. The long-range one- 
pion exchange transition potentials between isobar channels are important to this spectrum at low 
energies. We include them here in the 1 S'o channel and obtain good results for both elastic AW 
scattering and pion production. 

In previous work some isobar channels differing only in spin were combined. They are now 
treated as separate channels to more accurately treat the pion exchange potentials and for future 
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predictions of 3-body final state distributions. The effect of isobar width is included in more 
channels. 

In some instances the one-pion exchange (OPE) coefficients for particular channels were incor- 
rect in previous results. These corrections will be noted in the description of each partial wave, tq 
and 9 NN have been chosen so that one value is used throughout. 

In Section II we review critical aspects of the model. Section III describes the model elements 
and parameters for each NN partial wave and displays the fit to the phase parameters. In Section IV 
all the model partial waves with J < 4, and the one-pion-exchange amplitude for higher J partial 
waves are used to predict the experimental observables. The excellent results are displayed with 
full angular distribution of all spin observables for which there is data at key energies spanning 
Tl < 800 MeV. In addition, we present excitation curves at 90° for these observables and total 
cross sections (with and without polarization) over the whole range. In particular, we note that the 
well-known structures 9 in Aa L (pp) and Aa T (pp) are well- fitted with the exception of the somewhat 
shifted position of the Aa T peak. This defect is correlated with a slightly low minimum value of 
V ( 1 -^ 2 )- T ne physics that may account for this will be discussed. 

II. The Hadron Interaction Model 

A. The Schrodinger Equation 

For the internucleon distance r > ro, the coupled system wave function is determined by a 
homogeneous boundary condition at ro (to be discussed later) and a coupled channel Schrodinger 
potential with a meson exchange potential matrix. Using the strong interaction symmetries we 
assume conservation of total angular momentum, J, isospin I and parity. As is well-known this 
implies total spin, S, conservation in the NN system (neglecting the n — p mass difference). Con- 
sequently, the coupled system contains only one (L = J) or two (L — J ± 1, S — 1) NN partial 
waves, and may have as many NA, A A or AW* (1440) partial waves as have the same J, I and 
parity. Because of the width of the isobars each isobar channel partial wave is treated as many 
channels, 8 with the mass distribution discretized and weighted by the Breit-Wigner distribution. 
For each channel, i, the eigenvalue is determined by the relativistic value of the relative momentum 
ki in that channel, where the total energy W = y/Mf 2 + kf + yjAff + kf with Mf and M\ being 
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the pair of baryon masses in channel i. The system of Schrodinger equations is described in detail 
in Ref. [8]. 



B. The Potential in the NN Sector 

Within the NN sector, the potential matrix (diagonal and tensor coupling elements) is de- 
termined by the one-boson exchange of n, rj, p and uj mesons (Fig. 1(a)) and two-pion exchange 
(Fig. l(b)-(c)) in the adiabatic limit. These potentials, with the exception of the spin-orbit term, 
are described in Ref. [5] . In particular, the form of the two-pion potential is detailed in that reference 
including the "pair suppression" factor and the Breuckner-Watson vs. Taketani-Machida-Ohnuma 
ambiguity which arises in the adiabatic limit. All the coupling constants (with the exception noted 
below), the meson masses and the two parameters which arise in two-pion exchange potential are 
as in Ref. [5]. The exception is a small change in the value of the pseudoscalar coupling g^Nn- As 
a result of the larger core separation radius now used, and the coupling to isobars, the value now 
used is g% N7T /An = 14.40 (in place of the value 14.94 in Ref. [5]). 

The spin-orbit potential, L ■ S Vls, was ignored in Ref. [5] because of its short range. The tail of 
this potential beyond tq has been found to be of some importance to the L ^ partial waves. The 
sources of this term are two-pion, p and uj exchange, so that it has contributions of different range, 
form (the two-pion part is non- Yukawa) and sign. Because of the relatively small overall effect of 
this part of the potential, the radial form has been simplified by fitting a single Yukawa potential 
to the complete, relativistic result of Ref. [1 1] . In the 1=1 states this results in 

Vls = ~7 • (2.1) 

r 

It is smaller in the 1 = states and we neglect it. 

C. The Transition Potentials 

The off-diagonal components of the potential matrix connecting to isobar channels are obtained 
from one-pion exchange as shown in Fig. l(d)-(e). The general form is 12 

Vt{t) = ^ i • f 2 • S 2 v (r) + Sf 2 v 2 (r)} (2.2) 



where the fi are the coupling constants at vertex i, the Si (T^J are spin (isospin) or transition spin 
(isospin) operators at vertex i that will be defined below, the tensor operator is 



v (r) 



Sn = 3 (5i • f) (s 2 ■ r) - S! ■ S 2 , 

, ^ (, 3 3 
and v 2 (r) — 1 + 



m^r (m n rY 



v (r) 



(2.3) 
(2.4) 



The matrix elements of the above operators are expressible as 12 ' 13 
S'iS^S' L r J Si ■ S 2 SiS 2 SLJ^. 



(-1) 



S'i+S*2+5' 



Sss'S 



LL' 



S S' 2 S[ 
1 Si s 2 



l (spills 1 !) (^p 2 |5 2 ) 



(2.5) 



the analogous expression for T\ ■ T 2 and 



nt nf o f T ' T oT o O O T T 

Oio 2 o L J Oi 2 oio 2 oLJ 



^S"+J 



30(2L + 1)(2L / + l)(2S + 1){2S' + 1) 



1/2 



X 



J S' L'\ / L' 2 L\ I s ' 1 Sl 1 
2 L S | V 



)\s 2 S 2 1 |(^i||3i||5i)(^ 2 ||5 2 ||52) 
' { S' S 2 J 



(2.6) 



where the standard notation is used for the 3j (Wigner), 6j (Racah) and 9j coefficients. The 
Wigner-Eckart theorem requires 13 



S'rri 



S 



M 



S,m) = N (S'rri \\M,Sm 



(2.7) 



where the normalization N is determined by convention and the (S'm'\JM, Sm) are Clebsch- 
Gordan coefficients. In the Edmonds convention 14 N = 1 for the transition spin operators. For 
those cases diagonal in spin magnitude, because the Pauli spin operators and the equivalent spin- 
3/2 operator are twice the standard angular momentum operator normalization J = J{J + 1), 



N = -2 X /S(S + 1). For NNn, NN*n, iVAvr and AAvr vertices this requires 

-■Ve , 



2 IHI 2 



S 



and 



1/2,3/2 



S 



3/2,1/2 



= 2Vl5 



(2.8a) 
(2.8b) 

(2.8c) 



respectively, with analogous results for the isospin operators. 
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With these normalizations of the operators we have 

f NNn = 2M g NNn ( 2 '^ 

where in pp scattering, is the neutral pion mass and M is the proton mass. As discussed in 
Ref. [1], we obtain /nAty from the decay width of the A, giving ffj^/^ = 0.35 which is 50% larger 
than the quark model value but only 5% smaller than the result of a strong coupling model. 12 
Similarly we obtain 1 fNN*n from the iV*(1440) decay to a single pion, giving /| 0V » 7r /47r = 0.015. 
Lacking other information we have used the consistent quark and strong coupling model result 12 

5 (2.10) 



fl A j4n = 0.003 . 



D. The Diagonal NA Potential 



The equations and matrix elements for the diagonal NA potential (see Fig. l(f)-(g)) are of 
the same form as above in Section IIC. Here we discuss the implications of the fact that there 
are contributions from two distinct Feynman diagrams: the direct interactions, Fig. 1(f), and the 
exchange interaction, Fig. 1(g). To obtain its contribution to Vna,na, the exchange amplitude must 
be multiplied by (-1) L+S + T+1 . 

First we note that the coupling constant product for the direct interaction is fNNnfAAn/^ 71 — 

0.0155, while that for the exchange interaction is /nAty/^ = 0.351. However, that apparent 

dominance of the exchange interaction is largely reversed by the difference of the spin matrix 

elements. In the isotriplet state the matrix element of f ■ T is —5, but the matrix element of 
— » — » -1 

^1/2,3/2 ' ^3/2,1/2 is — 3- 111 addition, the spin-spin and tensor matrix elements are in most states 
larger for the direct interaction. The result is that the two contributions are of the same order, and 
as they sometimes differ in sign, the strength and sign of the long-range diagonal NA interaction is 
very sensitive to the channel involved. One may then hope to determine the effect of the diagonal 
coupling by comparing the effect in different channels. As we shall see, this is in fact the case 
in the 1 L>2 and channels, where strong inelasticity enhances the sensitivity to the diagonal 
components. 
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E. The .R-Matrix Boundary Condition 

As shown in Ref. [4], the boundary condition at tq is of the general form (a, (3 are channel 
indices) 

R^ = Y.MW)^{r) (2.11) 

where W is the total barycentric energy and the f a p(W) are meromorphic functions of W with poles 
of positive residue on the real axis. As discussed in Refs. [4] and [15], these poles are determined by 
the complete set of states in the interior satisfying ^i n t( r o) — 0- The pole position is given by the 
eigenvalue Wi of the corresponding interior state, and the residue by the derivatives of the interior 
wavefunctions at tq. It was shown in Ref. [16] that the residues can be re-expressed in terms of the 
rate of change of eigenvalue with r$ 

dWi ■ ■ 

P l afl = -n-a^S® ( 2 - 12 ) 
where the are the fractional parentage coefficients of the channel a in the interior state i. 
We then have 

AtfW=/& + E™^ ■ (2-13) 

i 1 

The fractional parentage coefficients relevant to our models are those of the (^i^) quark 
configuration which overlaps with the NN (^Soj, AA^So), NA^S^j and AA^S^j channels. 
They are £% N = |, £ AA = ^ A = ^ and ^ AA = -^=. 

References [2] and [3] discuss the theoretical constraints on ro, and the experimental determina- 
tion tq = 1.05 fm when the quark interior is described by the Cloudy Bag Model (CBM) parameters. 
In the CBM, in contrast to the MIT Bag Model, the experimentally determined ro is consistent 
with the theoretical constraints. 

III. Isobar Channels, Meson Exchange Potentials, Quark Configurations and Param- 
eters of the pp Partial Waves 

For J > 4 the amplitude used is that of OPE and the Coulomb force. Below we give the details 
of the model for each NN, 1 = 1 partial wave with J < 4. 

There are several properties of the model that are common to all the partial wave models. 

r / \-|6 

The i?-matrix separation radius, ro, is determined by requiring that the q ( IS1/2 ) configuration 
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energy, as determined by the CBM, be the energy for which the inner zero of the external wave 
function is at tq. The external wave function is determined by the potential matrix and the S- 
matrix. The data determines the latter for < 1 GeV. For the 1 S'o and s S\- 3 Di states the CBM 
and the data require 2 tq = 1.05 fm. The same value is indicated for the 1 L>2 state. 3 This value is 
about 20% smaller than the value of r corresponding to the equilibrium radius of these six-quark 
bags, which satisfies the i?-matrix method requirement that tq be within the region of asymptotic 
freedom. 2 We use tq = 1.05 fm for all partial waves. 

The coupling constants used at the meson-baryon vertices are those given in Refs. [1] and [5], 
with the exception that we use g^ N ^/A-K = 14.40 in all partial waves. This change is related to the 
larger tq now being used. One other consistent change from the earlier published work is an extra 
factor \f2 in the coefficients of the NN — > NA potentials. This is required by contribution of the 
NN — > AiV diagram, which is of the same strength as that of the NN NA diagram. 

Table I presents the values of all the isospin-isospin, spin-spin and tensor matrix elements that 
are needed for the determination of the OPE potentials between NN and isobar channels or between 
isobar channels, for isobar channels with one or two A's. The specific numerical coefficients are 
given in the subsection for each partial wave. 

We note that the /-matrix parameters determined by fitting the data depend slightly on the 
number and distribution of the discretized mass spectrum of the isobars. We have used 17 channels 
for each isobar whose width is not neglected. Other aspects of the distribution are discussed in 
Ref. [8]. 

A. The NN (^S ) System 

The last published versions of the coupled-channel i?-matrix model for this system were the 
"new S"' model of Ref. [3] and the "old S" model. 2 ' 3 Those versions were the first to include the 
CBM pole in the /-matrix with the associated larger value tq = 1.05 fm. The "new S n version 
takes isobar width into account. It was consistent with the experimental value of a [pp — > ppn + n~^j 
at Tl = 800 MeV. Its only important discrepancy from experiment was that it predicted too much 
inelasticity for 600 MeV < T L < 800 MeV. It predicted ??(800MeV) = 0.84 while the phase shift 
analysis (PSA) prefers 0.98. The excess of inelasticity was a result of the strength of coupling to the 
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Table I. 

Spin-spin and tensor matrix element for NN-NA, AW- A A, AA-AA and AA-AA transitions. 
a and a' represent the sets of quantum numbers in the first two columns. 
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(a'|5i-5 2 |a) 

\ | J. ^ | / 


(a/\SZ\a) 

\ 1 1Z \ I 
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3 
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+V6 






ii 












o 
o 

2 


a 
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2 
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2 


3 
2 


3 
2 


2 





V2 





2 


2 


3 
2 


2 








3 
2 


3 
2 








+V3 






ii 









2 


3 
2 


3 
2 


2 


-2V3 


V3 


1 





1 

2 


3 
2 


1 


1 





3 
2 


3 
2 


1 


-2 V / V3 







2 


, 

2 


3 
2 


2 





2 


' 

2 


3 
2 


2 


3 


-6 






ii 









2 


3 
2 


1 

2 


2 


+ 1 


1 





1 


1 

2 


1 
2 


1 





1 


1 
2 


3 
2 


1 


+ VW3 








ii 









1 


3 
2 


3 
2 


1 


—v/IO/3 


lv/275 





1 


1 

2 


3 
2 


1 





1 


1 

2 


3 
2 


1 


-5 


2 






ii 









1 


3 
2 


1 

2 


1 


+3 


5 
3 


1 


1 


, 

2 


2 


1 


1 


1 


2 


3 
2 


1 


+ VW3 








ii 






1 


1 


1 
2 


3 
2 


2 





3^3/10 






ii 






1 


3 


1 

2 


3 
2 


2 





3/V5 






ii 






1 


1 
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2 


3 
2 


1 


-v/TO/3 
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1 


1 


:. 

2 


3 
2 


1 


1 


1 


1 

2 


3 
2 


1 


-5 


-1 






ii 






1 


1 


3 
2 


1 

2 


1 


+ 1 


5 
6 


1 


1 


1 

2 


3 
2 


2 


1 


1 


1 

2 


3 
2 


2 


3 


21 
5 






ii 






1 


1 


3 
2 


1 

2 


2 


+ 1 


7 

10 


1 


3 


1 

2 


3 
2 


2 


1 


3 


1 

2 


3 
2 


2 


3 


24 
5 






ii 






1 


3 


3 
2 


1 

2 


2 


+ 1 


4 
5 


2 


1 


2 


: 

2 


1 


2 


1 


2 


3 
2 


1 


+ VW3 


+l/(5\/6) 






ii 






2 


1 


1 

2 


3 
2 


1 





-1^372 
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2 


3 


1 

2 


3 
2 


1 





3 
5 
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2 


3 


1 

2 


3 
2 


2 
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Q 
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1 

2 


1 

2 


1 
1 


Z 
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i 


1 

2 


3 
2 


1 
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5 
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2 


1 


1 
2 


3 
2 


2 





3 
5 
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2 


3 


1 

2 


3 
2 


1 
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2 


3 


1 

2 


3 
2 


2 





6 
5 


2 


1 


1 

2 


3 
2 


1 


2 


1 


1 
2 


3 
2 


1 


-5 


1 
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2 


1 


3 
2 


1 

2 


1 
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Table I. contd. 
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1 
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2 


1 
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3 
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2 


1 

2 


1 
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3 


1 
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3 
2 


2 
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3 
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2 
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2 


2 


3 
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2 


2 
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2 


1 

2 


2 
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5 


2 


2 


1 

2 
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2 
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2 
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2 


2 
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2 
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2 
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-2V3 






ii 






2 
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2 


3 
2 
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1 

2 
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2 


2 
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1 

2 


3 
2 


2 
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2 


1 

2 


2 
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2 


2 


1 

2 
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2 


2 


2 


2 
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2 
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3 
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3 


1 

2 


1 

2 


1 


3 


1 


1 

2 


3 
2 


2 
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3 
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3 


3 


1 

2 


3 
2 


1 
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3 


3 
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2 
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2 


2 





V3/10 


3 


1 


1 

2 
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2 


2 
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1 


1 

2 
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2 


2 


3 
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1 
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2 
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2 
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1 
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3 


3 


1 

2 
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2 


1 


3 


3 
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2 
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2 


1 


-5 
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3 


3 
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2 


1 

2 


2 
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5 


3 


3 
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2 


3 
2 


2 


3 


3 


1 
2 


3 
2 


2 


3 


14^2/9 






ii 






3 


3 


3 
2 


1 

2 


2 


1 


-7V2/27 


4 


3 


1 

2 


1 

2 


1 


4 


3 


1 
2 


3 
2 


2 
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4 


3 


1 

2 


3 
2 


1 




1/(3V6) 


4 


5 


1 

2 


1 

2 


1 


4 


3 


1 

2 


3 
2 


2 





-v/273 
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4 


3 


1 

2 


3 
2 


1 





-5^10/3 


4 


3 


1 

2 


3 
2 


2 


4 


3 


1 
2 


3 
2 


2 


3 


3 






ii 






4 


3 


3 
2 


1 

2 


1 


+ 1 


i 

2 


4 


3 


1 

2 


3 
2 


1 


4 


3 


1 
2 


3 
2 


1 


-5 


1 
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ii 






4 


3 


3 
2 


1 

2 


1 


+ 1/3 


5 

18 


4 


4 


1 

2 


1 

2 





4 


2 


1 
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2 


2 


3 
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7 


4 


2 


1 
2 


3 
2 


2 


4 


2 


1 
2 


3 
2 


2 


3 


12 
7 
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AA channel, which in turn was required to produce the correct energy dependence of the real phase 
5. Switching coupling strength to the high threshold AA or AW* channels reduced the inelasticity, 
but produced excessive two-pion production as in the "old S" model. 

The OPE potential in the diagonal AA interaction, now included, is repulsive, decreasing the 
inelasticity near the AA threshold. Because of the importance of two-pion production we also 
include the OPE transition potentials between the AA and AA states. We can now transfer 
AW — AA coupling strength to AW — A A coupling, but the increased production of these two-pion 
channels is transferred back to one-pion channels by the AA — AA coupling. 

The AA ( 5 A)) channel is now included along with the AA ( 1 5'o) channel, because of the 
former's strong OPE coupling. But the width of this channel is neglected, because its orbital 
angular momentum barrier means that its effective threshold is higher than its mass threshold. On 
the other hand, the AW* ( 1 >S'o) channel has been dropped because of its weak OPE coupling to the 
AW channel. It would be difficult to distinguish its effect from that of the AA (^Sb) channel. The 
result is a good fit to both 8 and 77 over the whole energy range. Simultaneously the production of the 
two-pion producing A A channel is kept low enough so that the contribution to o [pp — > pp7T + 7T~^J 
is about one-third of the experimental total value 10 of 3 fxb for Tl = 800 MeV. 

The potential matrix elements involving isobar channels are (Eq. (2.2)) 

V T [AW (^o) , AA ( 5 D )] = 0.31WU2 

= 0.172771 7r t>o 
= -0.1744771^2 

= —0.0696771^^2 

= 0.139771^770 — 0.069677i 7r t>2 and 
= -0.03837j + 0.1935tj 2 . 



V T [NN^So) ,AA 
AA( 5 A)),AA 



Vn 



V T [AA ( 5 A)) , AA ( 5 A) 

v D [aa ( 5 A) 



We note that the direct and exchange contributions to Vd (Figs. 3a and 3b) cancel in the 
coefficient of tjq, but add in the coefficient of V2 producing a substantial repulsive potential. 

The pole of the + , (ltSi/2) configuration is at W p = 2.71 GeV/c 2 , corresponding to 
T L = 2.04 GeV. Only the AA (^o) and AW (^o) channels overlap with the configuration, and 
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they have residues (Eq. (2.12)) 

Wiv = °- 149GeV > 

^iViV,AAO = °- 133GeV alld 

Paao,aao = °- 1 18718 GeV • 
Here and below we use the total spin S of the channel in the subscripts to distinguish otherwise 

ambiguous channels. 

The constant, non- vanishing, elements of the /-matrix are adjusted to fit the phase parameters 
and the singlet scattering length a p s p = — 7.82 fm. They are /$jvjvjv = 83.5743, /jva JVA = 0.0, 

/aao.AAO = 3 -°> /z\A2,AA2 = 5 -°> /jvjv,jva = 6 -° and /jvjv,AA2 = ~ 25 - 2 - The resulting phases 
are compared with the PSA results 17 ' 18 in Fig. 2, showing a very good fit. The PSA results of 
Refs. [19,20] generally agree with those of Refs. [17,18] but Ref. [19] does not go to the lower 
energies, and Ref. [20] has oscillations not present in the other PSA results or in ours. 

B. The NN ( 3 P ) System 

In the last published version 8 the NN ( 3 Po) channel was coupled to the NA ( 3 Po) channel only. 
We now also include the AA ( 3 Po) but neglect the effect of its width. In Ref. [8] the phase shift fit 
was improved for Tl < 300 MeV by reducing the pion coupling to g% Nlv /4ir = 13.0. However, with 
the increased value of tq used now (1.05 fm instead of 0.75 fm) we fit well with g 2 NN ^/An = 14.4 as 
in all the other partial waves. The added diagonal OPE potential for AA is repulsive. 
The potential matrix elements are 

V T [NN ( 3 P ) , AA ( 3 P )] = 0.207577^0 + 0.1037m^ 2 

V T [NN ( 3 P ) , AA ( 3 P )] = 0.130771^0 - 0.052777^2 and 

V D [NA ( 3 P )] = 0.1418777^0 + 0.0135777^2 . 

As seen in Fig. 3, a very good fit to the PSA is obtained with the non-zero /-matrix elements 

/ NN,NN = ll-?, / JVA, JVA = 2 - ' /aA,AA = 2 -°' /jVJV,JVA = 12 ' anci /jVJV,AA = 

We note that, both in our model results and the most recent PSA, the value of rj monotonically 
decreases below 1 GeV, unlike those of Ref. [8] . This is in part due to the repulsive diagonal OPE 
potential in AA. The result is that the previous indication of a resonance in the Argand plot (a 
counter-clockwise motion) is now gone for < 1 GeV. 



C. The NN ( 3 Pi) System 

Reference [1] contains the last published version of the NN ( 3 Pi) system. It uses the smaller 
tq and was coupled to a single AA(P) channel. The width of the NA channel could then only be 
taken into account when Vr was ig nored. We now couple to the NA ( 3 Pi), NA ( 5 Pi), NA ( 5 Pi) 
and the AA ( 3 Pi) channels. The width is taken into account only for the AA ( 5 Pi) channel, which 
is strongly core coupled as well as having a low angular momentum barrier and a low mass threshold. 
The other additions and changes common to all partial waves have also been made. The resulting 
potential matrix elements are 



V T [NN ( 3 Pi) , AA ( 3 Pi 

V T [NN ( 3 Pi) , AA ( 5 Pi 

V T [NN ( 3 Pi) , AA ( 5 Pi 
V T [NN ( 3 Pi) , AA ( 3 Pi 
V D [NA ( 3 Pi 
V D [NA ( 5 Pi 
V D [AA ( 5 Pi 



0.2075m^ - 0.05186771^2 + 0.8- 



,-2mir 



-0.2087777^2 + 3.0- 



r 

3 -2m,r 



= 0.1704777^2 - 2.5- 



0. 13777^^0 + 0.026777 7r 7;2 
0.1418777^0 - 0.00675777^2 
-0.1163777^0 + 0.0809777^2 
-0.1163777^0 + 0.924777^2 . 



In this case, fitting the PSA required the addition of the phenomenological two-pion range transition 
potentials to the OPE components connecting AW to AA channels. These phenomenological com- 
ponents substitute for p-exchange tails as well as actual two-pion exchange. They usually, as here, 
decrease the coupling strength of OPE. Although their numerical coefficients are large, because of 
the difference in the r-dependent form, they are dominated by the OPE tensor terms. 

The fitted non-zero /-matrix elements are (as before the numerals 1 and 2 in the subscripts 
denote the spin for P states, while the 3 denotes the F state) f^NNN = 17.915, /jvAl NAi = 2-0, 

fNA2,NA2 = 10 A /jVA3,iVA3 = 4 -°> /aA,AA = 5 -°> ZiViV,iVA2 = - 4 -°> and //V/V,AA = 10 - 85 - As shown 

in Fig. 4, the fits to PSA are very good. 
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D. The NN ( 3 P 2 - 3 F 2 ) System 

The last published version of the NN ( 3 P 2 — 3 P 2 ) system 8 already includes the effect of the 
width of the isobars. In addition to the general changes made to OPE couplings and the diagonal 
NA potential, the iVA 3 P 2 , 5 P2, 3 -^2 and 5 P 2 channels are each now treated independently, without 
combining the two P states and the two F states. The effect of width in the two F states is 
neglected because of the angular momentum barrier and the vanishing /-matrix coupling to the 
NN states. The fit to the PSA is better than in Ref. [8] and no phenomenological two-pion range 
potentials are now needed. The potentials involving isobars are 



V T [nn ( 3 P 2 
V t [NN ( 3 P 2 
V t [NN ( 3 P 2 
V t [NN ( 3 P 2 
V t [NN ( 3 P 2 
V t [NN ( 3 P 2 
V t [NN ( 3 P 2 
V t [NN ( 3 P 2 



, NA ( 3 P 2 
, NA ( 3 P 2 
, NA ( 5 P 2 
, NA ( 5 P 2 
, NA f 3 P 2 



, NA ( 3 P 2 
,iVA( 5 P 2 
,iVA( 5 P 2 

v D [na ( 3 p 2 
v D [na ( 5 p 2 

V d [NA ( 3 P 2 
V d [NA ( 5 P 2 



0.2075wo + 0.01037v 2 

-0.0762v 2 

-0.0933w 2 

0.0762v 2 

-0.0762v 2 

0.2075vo + 0.0415w 2 

0.1867t; 2 

0.1525w 2 

0.1418wo + 0.0014w 2 
-0.1163m - 0.0809w 2 
0.1418vo + 0.0054w 2 
-0.1163^0 + 0.02314w 2 



The /-matrix choice to fit the PSA is, for the non-vanishing elements, /jvtvp nnp = 13.917, 

/ NNF,NNF = 12 A fNAPl,NAPl = ^-5, InAP2,NAP2 = 2 -5> / NAF1,NAF1 = 4-0, / NAF2,NAF2 = 4 -0, 
f NNP,NNF = -2.8, / NNP,NAP1 = 10-30, and fMNP,NAP2 = — 1-2- 

The fit to the PSA shown in Fig. 5 is a substantial improvement to that of Ref. [8]. For two 
phase parameters it is not as good a fit as we achieve in most of the partial waves, but is adequate 
in most respects. There is little significance to the large difference between our 2 and that of Arndt 
et a/. 18 That parameter is very sensitive to the data and is strongly constrained in Ref. [18] by their 
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assumption of only one coupled isobar channel. 



graph configuration 



E. The NN ( 1 Z) 2 ) System 

In Ref. [3] the /-matrix pole arising from the 2 + state of the q ( 1 >S'i/2 / 
was included. The configuration energy, as given by the CBM quark dynamics and ro = 1.05 fm, is 
W p = 2.880 GeV. The configuration contains components of the AA (^S^j and AA (^S^j systems, 
and the non- vanishing residue-matrix components are (Eq. (2.12)) 

= 0.260 GeV 



AA2,AA2 



0.116276 GeV 



and 



^AA2,AA2 
PAA2,AA2 = °- 052 GeV • 

The rest of this previous model was described in Ref. [8], where the width of the A is taken into 
account. The OPE coupling to the A A {^S-^j channel was neglected then. 

Because its vanishing angular momentum barrier compensates for the higher threshold, we now 
include the AA ( 5 S 2 ) channel as well as the NA ( 5 D 2 ) and AA ( 5 S 2 ) channels included in Ref. [8]. 
The width of the AA channel is neglected. Among the general coupling changes made, the most 
important here is the inclusion of OPE diagonal potentials in the AA channels. They are attractive, 
greatly improving the shape of the inelasticity parameter, as will be shown below. The required 
potentials are 



V T 


~NN ( 1 D 2 ) 


NA ( 5 S 2 )~ 


= 0.1392^2 


V T 


NN ( 1 D 2 ) , 


NA ( 5 L> 2 ) 


= -0.1663v 2 


V t 


NN ( 1 D 2 ) , 


A A ( 5 S 2 )~ 


= -0.0794^2 




V D 


AA ( 5 S 2 )" 


= -0.0383vo and 




V d 


AA ( 5 £> 2 )" 


= -0.0383vo - 0.0415^2 



There are no phenomenological two-pion range potentials. The adjusted /-matrix elements are 

/ AW, AW = 15.60, /nAS,NAS = / NAD,NAD = 10 -°> /aA,AA = InA,NAS = ancl /Aw,AA = 
-6.50. 

The resulting fit, shown in Fig. 6, is quite good for 5. For 77, in contrast to the fit of Ref. [8], 
the rapid drop to Tl = 650 MeV followed by a sudden leveling is reproduced. The remaining defect 
is that the model value of rj in the 700 MeV region is too small by about 0.06 and the model value 



increases beyond 800 MeV. It is the attractive OPE diagonal NA potential which increases the 
effect of coupling at the lower energies, but is less effective at the higher energies, as required. 

F. The NN ( 3 F 3 ) System 

The Ref. [8] version of the NN (^F^j system, the last previously published, has the opposite 
deficiency to that just encountered in the 1 Z?2 state with respect to the inelasticity. The experimental 
r\ keeps decreasing uniformly for 515 MeV < Tl < 800 MeV, while the model rj begins to level off 
near 650 MeV. The OPE diagonal interaction in the NA (^Psj channel is dominantly repulsive (see 
below). The rapid drop of rj above threshold can be maintained by increasing the coupling to the 
NN channel, and the drop is then maintained to higher energy as the long-range repulsive potential 
becomes less effective. 

We now also include the NA ( 3 -Pb) and NA ( 5 -Pb) channels. The width of all NA channels is 
considered. The potentials are 

' 1 = -0.1822^2 

= 0.2075wo - 0.0520^2 
= 0.0696^2 

= -0.1163uo + 0.0231^2 
= 0.1418wo - 0.0067^2 and 
= -0.1163^0 - 0.0424^2 . 
No phenomenological two-pion potentials are required. The non-zero fitted /-matrix parameters 

are fNN,NN = 5 - 90 > ZiVA2,iVA2 = °- 8 > J'nA1,NA1 = L0 > and J'nN,NA2 = As shown in Fi S- 7 - a 

very good fit is now obtained. 

G. The NN ( 3 F 4 - 3 if 4 ) System 

In the case of the NN ( 3 i<4 — ^H^j system the last published version dates back to the no-isobar 
results of Ref. [1]. In Ref. [3] it was remarked that coupling to isobars would not have a large effect 
for Tl < 1 GeV because of the large angular momentum barrier (L > 3). In fact, the change due 
to the couplings described below is not large, but there is a substantial improvement to the fit to 
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V T [NN ( 3 F 3 ) , NA ( 5 P 3 
V T [NN ( 3 F 3 ) , NA ( 3 F 3 
V T [NN ( 3 F 3 ) , NA ( 5 F 3 
V D [NA ( 5 P 3 
V d [NA ( 3 F 3 
V D \NA f 5 F 3 



The couplings to both the NA (^F^j and NA ( 5 -^i) channels are included, and the width of the 
A is not neglected in either channel. The potential matrix elements are 

V T [NN ( 3 F 4 ) , NA ( 5 F 4 )1 = -0.1 16^2 



V, 



0. 1032u 2 -2.1 e~ m * r /r 
0.207wo + 0.017^2 
-0.078^2 + 1.6 e- m * r /r 
-0.116^0 - 0.0578^2 and 
0.142vo + 0.0025^2 . 



T [NN ( 3 # 4 ) , NA ( 5 F 4 ^ 
V T [NN ( 3 F 4 ) , NA ( 3 F 4 } 
V T [NN ( 3 # 4 ) , NA ( 3 F 4 ^ 
V D [NA ( 5 F 4 } 
V D [NA ( 3 F 4 } 

The strongly attractive OPE diagonal potential in the NA (^F^j state compensates substantially 
for the effect of the orbital angular momentum barrier, to the extent that too much low-energy 
coupling occurs. Consequently, a best fit requires the inclusion of a phenomenological two-pion 
range contribution in the NN {^H^ coupling to NA (^F^j and NA {^F^ , as above. 

The fitted non-zero /-matrix constants are ZjviVF,iViVF = 8 - 49 > fNNH,NNH = 30 -°> /jVA2,JVA2 = 5 -°' 
ZiVAi,iVAi = 5 -°> f nnf,nnh = - 2 -°> f NNF,NA2 = -3.90, and Jnnf,nai = 4 - 80 - The results, shown 
in Fig. 8, are in good agreement with the PSA except for 5 (^H^j. The predicted values of 5 (^Ft^J 
are too attractive over most of the energy range. The coupling and the boundary condition have 
little effect because of the L = 5 barrier. If the PSA results are correct this appears to represent 
a deficiency of the potential in the NN sector for large L, perhaps indicating a repulsive quadratic 
spin-orbit term. 

H. The NN ( X G 4 ) System 

As in the previous NN ( 3 F 4 — 3 iJ 4 ) case, modified versions of the NN ( 1 C 4 ^) system have not 
been published since the uncoupled result of Ref. [1] . In this case the coupling is potentially more 
important because it couples to the L = 2 NA {^F)^j channel. The width of that channel is included. 
The isobar related potential matrix elements are 

V T [NN ( j G 4 ) , NA ( 5 D 4 )] = 0.1663^2 - 1.80 e^Vr and 
V D \NA ^dA] = -0.0383wo + 0.0553^2 . 



A phenomenological two-pion range contribution in the coupling potential is again required to cancel 
a coupling effect which is too strong at medium-range. The /-matrix constants are f^fNNN = 9-80, 
/jva AA = 2-5 and /jvtvaa = 4-0. As shown in Fig. 9 a very good fit to the PSA is obtained. 

IV. Amplitudes and Observables 
A. The Method 

With respect to the observables, several conventions and notations are in use in the literature 
and it is sometimes difficult to determine their relationship. This is only partly due to notation 
in which zero-spin labels were omitted and triple-scattering experiments were not considered or 
clearly distinguished. It mainly owes its origin to the different definitions of the coordinate systems, 
both in the center of mass (cm.) and in the laboratory. Here, we stick to the Madison convention 
(positive normal along ki x kf) and the Saclay convention 21 (same normal used in both frames but 
also in the three different sets of basis vectors in the laboratory) and avoid the Basle convention 
which uses two normals. Five common choices of orthonormal basis vectors are compared in Table 
II. 

Table II 



Group 


Center of mass system 


Laboratory system 


Saclay 21 


ki 




n 


£ 


m 





k 


s 


k! 


s' 


k" 


s" 


n 


Wolfenstein 22 


V 


P' 


n 


P 


K 


omitted 


k 




k! 


s 


K 


St 


n 


Halzen- Thomas 23 


z 




y 


t 


t 


omitted 


I 


s 


1 


s 


i 


s 


n 


Argonne 


same as above 


L 


S 


L 


S 


L 


S 


N 


Raynal 25 


ki 




n 


P 


K 


omitted 


z 


—x 










y 



A selection of some representative sets of basis vectors encountered in the literature both in the 
center-of-mass and laboratory systems. The vectors listed in a column are all equal. One notes that 
(f) £ = x sin | + z cos | and {\) m = x cos | — z sin | . 



The notation (Beam, Target; Scattered, Recoil) to express observables is adopted by Argonne 
whereas Saclay used (Scattered, Recoil; Beam, Target). The position of a given spin label in 
such a quadruplet thus specifies to which particle the frame is attached. Taking into account this 
information, the use of superscripts in the lab system by Saclay, in order to distinguish the incident 



(no superscript), scattered (prime) and recoil (double dash) particle is unnecessary and is being 
progressively abandoned in observable names. It follows from Table II that Ckk = Cmmoo = Aoo mm , 
Cpp = Cuqq = Aqqu and Ckp = Cpx = QmOO — ~~ ^ooim (i n eac h case the last equality comes 
from time reversal invariance). The Wolfenstein system 22 is then simply connected to the Saclay 
system. 21 

The theoretical group at Argonne is used to the Halzen-Thomas notation. 23 Their vectors 
(x, y, z) define the center-of-mass helicity frame. The experimental group at Argonne has capitalized 
the spin labels in the laboratory system. Raynal's notation 25 defines the spin labels x, y, z in the 
laboratory system but these should not be confused with the Argonne center-of-mass same labels x, 
y, z. The polarization labels X, Y, Z used in the SAID system correspond to Raynal's index x, y, 
z. Convention also differs, as the sideways direction x is chosen to be parallel to —s. Such a choice, 
however, is not consistent with the Saclay right-handed frame attached to the incident particle in 

— * 

the lab system and for which n — k x s. This is the reason why some minus signs appear when 
connecting Raynal's formalism 25 to the Saclay one, e.g. Azx = — A)0sfc- Still another different 
convention but now in the center of mass can be found in the literature. In some older papers, 
the defined vector m has its sign opposite to the one of the Saclay formalism. This was made so 
that in the non-relativistic case the cm. vector m coincides in direction with the lab momentum 
k" of the recoil particle, which is not the case in the Saclay formalism. Again connecting identical 
observables defined with respect to two different triplets of basis vectors in the same coordinate 
system involves a few minus signs (as many as there are in index). A final source of confusion is 
that the Wolfenstein transfer parameters R' t and A' t are defined with respect to —k", which changes 
the sign of the normal in the lab frame attached to the recoil particle. The connection with the 
Saclay formalism is then R' t = — i^ofe"sO an d A' t = —K^ii^. 

Keeping these five choices in mind, we compare the observables that we will compute, in Table 
III. One can see that the names used in the SAID system come from different groups. 

We calculate the observables from the bilinear combinations of the helicity amplitudes, which 
themselves are constructed from the partial-wave S'-matrix. The spin partial-wave expansions of 
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Table III 



SAID 18 


Wolfenstein 22 


Saclay 21 


Argonne 24 


DSG 


h 




a 


P 


P 


-PnOOO = ^OOnO = ^OOOn 


P 


D 


D 


DnOnO = DonOn 


Dnn 


DT 


D t 


KnOOn = KQ nn Q 


Knn 


R 


R 


Ds'OsO = Dq s "q s 


Dss 


RP 


R' 


Dk'OsO = D ok » 0s 


D S L 


A 


A 


D s '0k0 = A)s"0fc 


Dls 


AP 


A' 


Dk'OkO = A)fc"0fc 


Dll 


AYY 




CnnOO = Aoo nn 


Cnn 


AZZ 


Cpp cos 2 | + Ckk sin | 
+ Ckp sin 




Cll 


AXX 


Cpp sin 2 1 + C K k cos 2 § 
— Ckp sin 9 


AqOss 


Css 


AZX 


\{C K k - Cpp) sin6> 
+ Ckp cos 9 


-A()O s k = — A)0fcs 


-CSL 



Name and sign conventions for the observables computed in the present paper in 
relation to four different formalisms often encountered in the literature. 

the s-channel helicity amplitudes follow: 

01 = 2^ E { (2J + 1)a J,o + ( J + V a J,+ + Ja J,- + 2^(^ + 1)^} ^j(cos^) (4.1a) 

02 = ^ E {-( 2J + ^"j.o + ( J + + Jc V + 2 V /j ( J + 1 )« J } P J( C0S ^) ( 4 - 1& ) 

03 = ^ E {( 2J + l ) a j,i + Ja j,+ + ( J + l ) a j,- ~ 2\/j(J+lW} x 



x 



(4.1c) 



Pj (cos 9) + \. T COS _ i d . P! J (sm9) 



j(j + iy 



■4 = ^ E {-( 2J + ^"j,! + + ( J + ^"j- " 2^7(7 +l)a J j 



x 



(4. Id) 



x 



-^(cos^) + i±^P^(sin^) 



\ 00 c 



j| sin#P^(sin#) 



(4.1e) 



-20- 



where Pj{cos9) are the Legendre polynomials of order J (total angular momentum), Pj{sin9) 
are derivatives of Pj(cos9) with respect to their argument, and k is the cm. relative momentum. 
The energy-dependent a parameters describe the different possible S- and L- coupled states: a JQ 
represents the S singlet, a J1 represents the uncoupled triplet, and ol j± , a J represent the L coupled 
triplet states. It is understood in the above sums that, when we encounter a J = pole, then the 
corresponding a parameter has to be set equal to zero. Consequently, a Q 1 = a Q + = a Q _ = a = 
and the 3 Po lacks a coupled partner. The a J± are the diagonal partial- wave amplitudes for tensor 
coupled states while the a J are the off-diagonal ones. 

If S is the real part of the complex phase shift and ^/fj is the exponential of the imaginary part, 
then our inelasticity parameter 77 is related to Bugg's p [Ref. 17] by the equality rj = cos 2 p. Thus 
77 = 1 in the elastic region Ti < 300 MeV. For the parametrization of the partial-wave amplitudes 
in terms of the phase shift parameters, that definition leads to 

a j,o = ^j,o ex v( 2i6 Jfi) - 1 ( 4 - 2a ) 

= V ex P( 2 ^J,i) - 1 ( 4 - 26 ) 
a J± = rj J± cos2e J exp(2i5j ; ±) — 1 (4-2c) 

qJ = i \f*h^\+ sin 2I J exp [* + 5j >+ + ( 4 - 2rf ) 

where r] J± = 1 in the elastic region and €j is the elastic (or "real") mixing angle between J, + and 
J, — states. The inelastic (or "imaginary" ) mixing angle (pj is non zero only above the inelastic 
threshold and its numerical value remains small just above it. However, for states of small angular 
momentum barrier {e.g. 2 ), its value rises rapidly with energy and has a point of inflection at 
the A resonance threshold. The mixing angles 2 and 4 presented in relation (4. 2d) should not 
be confused with the helicity amplitudes of eqs. (4.1b) and (4. Id). In that representation, the e's, 
5's and 0's are real numbers and < 77 < 1. In our figures for 77's and <5's, we replace the set of 
label subscripts by the spectroscopic notation { 2S+1 Lj). The coupled channel parameters obey a 
sub-unitarity inequality 26 because of their other coupling to isobar channels. Finally, the 3 Po wave 
has the same representation as the singlet states. 

The observables can be computed from the cm. invariant amplitudes used by Saclay (formulae 



can be found in Ref. [21]) and related to the helicity amplitudes by 

a = \ [(^! + 02 + 03 - 04) cos 9 - 40 5 sin 9] (4.3a) 

b = l - (0i - 02 + 03 + 04) (4.36) 

c= l - (-01 + 02 + 03 + 04) (4.3c) 

d = \ (01 + 02 - 03 + 04) (4.3d) 

e= ~2 f^ 1 + ^2 + 03 - 04) sin 9 + 405 cos 9} (4.3e) 

This particular set of amplitudes is transformed into several other ones also used in the literature, 
in Refs. [21] and [27]. The normalization is such that 

^ = ^{H 2 + N 2 + N 2 + H 2 + N 2 } (4.4) 

is in mb/sr. For J < 4, the partial waves used are the coupled channel ones of Section III, that is 
^o, 3 P , 3 -Pi, 3 ^2 - 3 F 2 , 1 J D 2 , 3 ^3, 3 F 4 - 3 H 4 and ^ In fitting the energy dependence of the total 
cross sections, a compromise is reached in using the waves and x Iq given by the Feshbach-Lomon 
model without isobar coupling. 4 The same NN coupling constant in the ir- and 2ir- sector (14.40) is 
used throughout. The unitarized OPE phases begin with 3 Hq — 3 Jq, 3 Jt, . . . until J = 12, the starting 
value for the OPE Born amplitude. The Coulomb interaction is included, except for the total cross 
sections obtained by the optical theorem. In the latter case, the pure Coulomb amplitudes and 
phases are removed but the electro-nuclear interferences in the model phases are kept. 

The above prescriptions for treating OPE and Coulomb effects in our model are of course slightly 
different than those of Bugg 17 or Arndt. 18 Nevertheless, we will see that the different treatments 
lead to minor differences in the polarization observables at low energy (142 MeV) and that the 
differences become appreciable only when 9 < 5° and only at higher energies. The observables of 
Table III are calculated at energies of 142, 210, 325, 400, 475, 515, 580, 615, 685, 750, 800 and 970 
MeV. The excitation functions to be displayed below cover the whole range of energies, Tl < 1 
GeV, whereas angular distributions are presented at 142, 515 and 800 MeV only (except for A zz 
for which we add 580 and 685 MeV). 
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B. The Results 

The excitation curves for the total unpolarized (er r ) and polarized (Aer^ and Aa" T ) cross 
sections are shown in Figure 10. The experimental points are in general well reproduced by the 
model predictions as the latter follow the correct curvatures and structures in the energy dependence, 
especially for Aa L . However Aa T does deviate from the data in that the position of the model peak 
is about 50 MeV (laboratory energy) higher that the data peak, and the prediction falls off more 
slowly that the data after the peak. This is correlated with the deviation, already noted in Section 
III, of i]( 1 D 2 ) from the PSA result for T L > 600 MeV. The model of i]( 1 D 2 ) continues to decrease 
after the PSA result levels off at 600 MeV, and it does not rise to the PSA value until Tl > 800 
MeV. Because Aa T increases with decreasing r]( 1 D2), the model peak is shifted to higher energy 
and the value is higher than the experimental value until 900 MeV. 

The excitation functions at 90° cm. for R, R' , A, A' , D, C nn oo, ^oofcfc; A)0ss and Ckp are 
compared with the data points between 89° and 91° compiled in the SAID system, 18 and are dis- 
played in Figure 11. Our predictions for Ckk are also compared with Arndt PSA results. For iden- 
tical particles at 9 cm = |, other theoretical predictions are P = 0, Aoosfc = and D n Q n Q = Ko nn o 
(the few existing Dt experimental points are thus plotted together within the D experimental set). 
We note that R' and A' are not independent quantities at this angle. For pp scattering at 90°, one 

A 7?' 

has -jj = — = tan#i with 9\ being the lab scattering angle. Theory implies that the polarization- 
transfer parameters here are not independent from the depolarization and rotation parameters since 
R t (90°) = R(90°), R' t (90°) = R\90°), A t (90°) = -A(90°) and 4(90°) = -A' (90°). Consequently 
we could have plotted any independent measurements of Rt, R' t , A t , A' t together with the usual 
Wolfenstein set, as we did for D t and D. However, at this time no data points exist for the transfer 
parameters in the range 88° < 9 < 92° up to 1 GeV. Under the same circumstances, one also has 
^oofcfc + Aoo ss = C nn oo — 1 because of the Pauli principle. Other constraints at 90° cm. on the 
limits of some spin observables from amplitude analysis can be found in Ref. [28]. 

The fit of the 90° excitation functions to the data is generally very good. The only substantial 
deviation from data for T L < 800 MeV is for ^(90°) at 650 MeV < T L < 800 MeV. In this range 
the newest data is higher than the older data and our prediction. Above 800 MeV, up to 1 GeV, 
our predictions for D, A yy and A zz decrease with energy while the data remains nearly constant. 
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For the El ', A', A xx and excitation functions we also show the result of the SAID PSA 18 
because the data is sparse or non-existent. The trend is similar to our prediction, with somewhat 
less structure. 

The comparison with the angular distributions is shown in Fig. 12. The model energies chosen 
as representative are T L = 142, 515 and 800 MeV for P, R, R', A, A', D, D t , A yy , A zz , A zx , 
and A xx . In addition we compare with recent A zz data (LAMPF) at 580 and 685 MeV. Again the 
Arndt PSA prediction 18 is also shown where the data is sparse or non-existent. The experimental 
energies compared with these curves are listed in the caption. They are mostly within 5 MeV of 
the model energy, but sometimes differ by as much as 10 MeV. There also exists some Wolfenstein 
spin-transfer R t and A t data, 18 but only at 800 MeV. There are eight points of moderate to poor 
accuracy. None of the points are near 9 cm = 90°. We have not calculated the model predictions for 
Rt and At or any of the triple-index spin observables. 

The fit to the angular distributions is good, and has no important deviations from the data for 
Tl < 800 MeV. Most of the observables are also well reproduced at Tl = 800 MeV. However at 
this energy the model shows more angular dependence than the data for the observables P, R', A, 
Ayy, and A zz . This trend persists at 580 and 685 MeV for A zz , but the statistical significance is 
small. The problem may be associated with the model value of 5( 3 i?4). This has the only significant 
deviation from the PSA of any of the 5's. The high L value implies a strong angular dependence, 
and our 5( 3 i?4) is significantly larger than the PSA for Tl > 500 MeV. 

V. Conclusions 

We have presented the results of a model for pp scattering that includes: 

(a) Coupled NN, NA and AA channels for J < 4. The effect of the width of the A is taken into 
account. 

(b) The 7i, 2n, rj, p and uj meson exchange potentials in the NN sector, and 7r-exchange potentials 
connecting the isobar channel sector to itself and to the NN sector. The coupling parameters 
used were determined by experiments other than NN scattering or from SU (3) and quark model 
relations. 
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(c) At a separation radius r$ an i?-matrix boundary condition is imposed representing valence quark 
degrees of freedom for r < r$. The energy independent part of the /-matrix are parameters. 
The poles are determined by the energies of the quark configurations as given by the CBM. 
Only the poles of the lowest 6-quark configuration (lsi) 6 are included, affecting the ^Sq and 
1 D2 PP channels. The value of tq is determined by the physical requirements of the i?-matrix 
method. 

The result is a very good fit to almost all the data for Tl < 1 GeV. The major exceptions are the 
higher energy of the peak in Aa T and an excessive angular dependence in a few of the observables. 
No other model produces quality fits for Tl > 300 MeV and they especially fail in describing the 
structures in Aa L and Aa T . The good degree to which these structures are fitted by this model 
strongly implies that the structures are due to the effect of coupling to the NA threshold in the 
NN 1 L>2 and 3 F3 channels. 

The shift of the Aa T peak is related to the fact that the model r)( 1 D2) structure is deeper than 
the PSA result. The inability of the model to more closely fit the PSA value of rj( 1 D2) may be 
related to the omission in the model of coupling to the 7rD channel. This channel couples relatively 
strongly to the NN( 1 D2). Including it, in addition to the isobar channel coupling of the model, 
may allow a simultaneous fit of the steep descent of rj( 1 D2) just above threshold and the weakened 
inelastic effect above 650 MeV laboratory energy. 

The excess of angular structure in some of the observables for Tl > 515 MeV is likely due to the 
model value of 5( 3 i?4) being larger than the PSA value at the higher energies. The difference is at 
most 0.5°, but this can be significant because of the high statistical weight of the J = 4 partial wave. 
Because of the very large angular momentum barrier (L = 5) the boundary condition parameter 
has negligible effect. This indicates that the medium range meson-exchange potential used here 
is inaccurate for large L. A correction to our potential can arise from high order non-local effects 
which give rise to quadratic spin-orbit type terms. These were not included. 

In spite of the above noted deficiencies the pp data for < 1 GeV is overall well understood 
in terms of this model implying that most of the physics has been included. An important part of 
this physics is the i?-matrix boundary condition representing the quark degrees of freedom at short 
range. Even in those channels in which we ignore the pole terms (because they are distant), the 



energy independent boundary condition acts differently than any energy-independent non-singular 
potential. 4 In fact it behaves as if the interior potential increases linearly with energy, asymptotically 
becoming a hard core 29 , and the present results reinforce that conclusion. 

Of course the most interesting effects of the i?-matrix core representation come from the poles 
coinciding with quark configuration energies. These give rise to exotic quark resonances, exotic 
dibaryons in the present case, which begin at Tl > 1.8 GeV. The implications for structure in 
the observables have already been given for the earlier forms of these model 2,3 and have interesting 
correlations with experiment. 3,30 ' 31 Those predictions will soon be updated with the present models, 
and the models themselves will be extended to include the poles due to the (lsi ) 5 (lpi ) configuration. 
This modification to the odd-parity states will negligibly affect the Tl < 1 GeV predictions. 

The other obvious extension of this work is to update the 1 = partial wave models and 
compare with the np data. That work is well underway. 
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FIGURE CAPTIONS 

Fig. 1: Meson-exchange potential diagrams, (a)-(c) Contributions to NN-NN sector potentials; (d) for 
AA-A-isobar potentials; (e) for NN-AA potentials; (f) for NA-NA potentials; (g) for NA-AN 
potentials. 

Fig. 2: Phase parameters for AA ( 1 5'o) scattering. Model - solid curve; PSA of Ref. [18] - dashed 
curve; PSA of Ref. [17] - □. 

Fig. 3: Phase parameters for AA ( 3 -Po) scattering. Notation as in Fig. 2. 

Fig. 4: Phase parameters for AW ( 3 -Pi) scattering. Notation as in Fig. 2. 

Fig. 5: Phase parameters for the tensor coupled AA ( 3 P2 — 3 ^2) scattering. The parameters determine 
the ^-matrix 

S— = V (V " 1) j) cos 2ej exp 2i5 ( 3 ( J - 1) j) , 

S+- =S-+ = i [v ( 3 ( J - 1) j) V (V + 1) j)] sin 2ej 

x exp [5 ( 3 (J- 1) J )+5( 3 (J + 1) J )+0 J ] , 

S++ = V (V + 1) j) cos 2ej exp 2i5 ( 3 ( J + 1) j) . 
Notation as in Fig. 2. 

Fig. 6: Phase parameters for AA i^D^ scattering. Notation as in Fig. 2. 

Fig. 7: Phase parameters for AA ( 3 -^s) scattering. Notation as in Fig. 2. 

Fig. 8: Phase parameters for the tensor coupled AA {^F/± — ^H^j scattering. Parameters determine 
^-matrix as in Fig. 5. Notations as in Fig. 2. 

Fig. 9: Phase parameters for AA i^G^j scattering. Notation as in Fig. 2. 

Fig. 10: Excitation functions for the total unpolarized (°" TOr ) and polarized (Aa L and Aa T ) pp cross 
sections for Ti < 1 GeV. 
□ : experimental values compiled in SAID. 18 

for Aa T : Saclay points (PE86). 
Solid curves: model. 
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Fig. 11: Excitation functions at 6 cm = 90° for the differential cross section Wolfenstein parameters 
R, R' , A, A' , D, second-rank asymmetry spin-tensor A yy , A zz , A xx (A zx = at this angle), and 
spin-correlations Ckp and Ckk- The two last parameters are defined in the cm. frame but 
the experimental points shown do represent independent measurements. Note that Cpp (90°) = 
G^(90°). 

□ : experimental values compiled in SAID, 
x for ^: Saclay points (GA85). 

x for D: independent D t points used since D(90°) = D t (90°) (they are listed in SAID 18 ). 
Solid curves: model. 

Dashed curves: Arndt's PSA predictions. 18 

Fig. 12: Angular distributions computed at = 142, 515, and 800 MeV compared with experimental 
values compiled in SAID 18 around each energy, for the following observables. 

^ at 137 (0), 144 (O) and 144.1 (x) MeV 

513 (□) and 516 (x) MeV 

788.7 (0), 789 (x), 795 (□) and 800 (O) MeV. 
P at 137 (0), 140.7 (□), 142 (O) and 147 (x) MeV 

515 (□), 515.3 (0) and 517 (x) MeV 

790.1 (0), 794 (x), 796 (□) and 800.9 (O) MeV. 
R at 140 (□) and 142 (O) MeV 

515.3 (□) and 517 (x) MeV 

800 (O,D,0) MeV. 
R' at 137.5 (□) and 140.4 (O) MeV 

515.3 (□) and 520 (x) MeV 

800 (x,D,0) MeV. 
A at 139 (□) and 143 (O) MeV 

517 (□) MeV 

800 (O,D,0) MeV. 
A' at 800 (x,D,0) MeV. 

-30- 



D at 138 (□), 142 (O) and 143 (x) MeV 
515.3 (□) and 517 (x) MeV 
800 (O,x,D,0) MeV. 
D t at 517 (□) MeV 

800 (0,D) MeV. 
A xx at 514 (□) MeV 
791 (O) MeV. 
-A zx at 514 (□) MeV 

793 (x), 794 (O) and 805.7 (□) MeV. 
A zz at 514 (□) and 518.4 (x) MeV 

577 (x), 583 (O), 586.3 (0) and 589 (□) MeV 
688 (x) and 692 (□) MeV 

790.1 (0), 793 (x), 800 (O) and 805.7 (□) MeV. 
A yy at 143 (□) MeV 
515 (□) MeV 

790.1 (0), 794 (x), 796 (□) and 800 (O) MeV. 
Curves notated as in Fig. 11. 
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For a copy of figure 1 

[the only one which is not in PostScript] 

please send an email request 

containing your mailing address 

to ereidell@marie.mit.edu 

Thank you! 

— Evan A. Reidell (MIT CTP support staff) 
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